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ABSTRACT 

We calculate non-axisymmetric oscillations of neutron stars magnetized by purely 
poloidal magnetic fields. We use polytropes of index n = 1 and 1.5 as a background 
model, where we ignore the equilibrium deformation due to the magnetic field. Since 
separation of variables is not possible for the oscillation of magnetized stars, we em¬ 
ploy finite series expansions for the perturbations using spherical harmonic functions. 
Solving the oscillation equations as the boundary and eigenvalue problem, we find two 
kinds of discrete magnetic modes, that is, stable (oscillatory) magnetic modes and un¬ 
stable (monotonically growing) magnetic modes. For isentropic models, the frequency 
or the growth rate of the magnetic modes is exactly proportional to Bg, the strength 
of the field at the surface. The oscillation frequency and the growth rate are affected 
by the buoyant force in the interior, and the stable stratification tends to stabilize the 
unstable magnetic modes. 
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1 INTRODUCTION 

Neutron stars are believed to have strong magnetic fields. The field strength Bs at the surface is estimated to be ~ 10^^ G for 
pulsars and ~ 10^® G for magnetar candidates, but we do not have any good knowledge of the configuration of magnetic fields 
inside the star. Recently, quasi-periodic oscillations (QPOs) in the tail of the giant X/y-ray flare were detected from the soft 7 
ray repeaters (SGRs) (e.g., Israel et al. 2005; Strohmayer & Watts 2005). SGRs belong to what we call magnetars. Giant flares 
have so far been observed only from three SGRs, that is, SGR 0526-66, 1900-1-14, and 1806-20, and just once for each of the 
SGRs, indicating that giant flares in magnetars are quite rare events. The QPOs are now regarded as a manifestation of global 
oscillations of the underlying neutron stars, and it is expected that they can be used for seismological studies of the magnetars. 
However, it is not an easy task to properly carry out modal analysis of magnetar candidates because of their extremely strong 
magnetic fields, which can significantly modify the modal property of the stars. Because separation of variables between the 
radial and angular coordinates is not possible for the perturbations in magnetized stars, we have to employ series expansions 
to represent the perturbations in linear analysis, which possibly makes the analysis difficult, particularly when singular points 
inherent to the governing equations appear in the interior of the star (see, e.g., Asai, Lee, & Yoshida 2015). 

Many authors have calculated axisymmetric (m = 0) oscillation modes of magnetized stars having a purely poloidal 
magnetic field, where m is the azimuthal wave number of modes. Using a toy model, Glampedakis, Samuelsson, and Andersson 
(2006) and Levin (2006, 2007) investigated global oscillation modes residing in the fluid core and in the solid crust both 
threaded by a magnetic field, and they showed that frequency resonance between modes in the core and in the crust could be 
important. Lee (2008) and Asai & Lee (2014) carried out normal mode calculations of axisymmetric toroidal modes and found 
discrete modes, van Hoven & Levin (2011, 2012) suggested the existence of discrete modes in the gaps between frequency 
continua, using spectral method. Besides the normal mode analyses mentioned above, however, most authors have used MHD 
simulations to investigate the modal properties of magnetized stars (e.g., Cerda-Duran et al. 2009; Golaiuda & Kokkotas 2011; 
Gabler et al. 2011, 2012; Sotani et al. 2008). Sotani & Kokkotas (2009) calculated axisymmetric polar-Alfven oscillations and 
found that continuous frequency spectra are not formed. Golaiuda & Kokkotas (2012) calculated mixed polar and axial 
oscillation modes of magnetized star having both poloidal and toroidal magnetic fields for axisymmetric perturbations. They 
suggested that the oscillation spectra can be significantly modified by the toroidal magnetic field component. Gabler et al. 
(2013a) calculated magnetic oscillation modes for various magnetic field configurations (e.g., purely poloidal, purely toroidal, 
mixed poloidal and toroidal). 

For non-axisymmtric (m 7 ^ 0) modes of magnetized stars, however, there are only a few numerical studies of global 
oscillations. For purely toroidal magnetic field configuration. Lander et al. (2010), Passamonti & Lander (2013), and Asai, 
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Lee, and Yoshida (2015) calculated non-axisymmetric oscillation modes. Lander & Jones (2011) calculated non-axisymmetric 
oscillation modes of magnetized rotating star with purely poloidal magnetic field using MHD simulations, and they obtained 
polar-led Alfven modes, which reduce to inertial modes in the limit of M/T —>■ 0, where M and T are magnetic and rotation 
energies of the star. It is important to note that Lander & Jones (2011) suggested that the axial-led Alfven modes could be 
unstable. 

In this paper, employing the normal mode analysis, we calculate non-axisymmetric (m 7^ 0) oscillation modes of magne¬ 
tized stars having purely poloidal magnetic fields. We assume that the gravitational energy dominates the magnetic energy so 
that the stellar deformation due to the magnetic field can be safely neglected. No effect of rotation is considered. For normal 
mode analysis in this paper, we employ finite series expansions for perturbed quantities to derive a finite set of coupled linear 
ordinary differential equations, which is solved as the boundary and eigenvalue problem by imposing appropriate boundary 
conditions. This paper is organized as follows. §2 describes the method used to construct a magnetized equilibrium stellar 
model, and perturbation equations for non-axisymmetric oscillation modes in magnetized stars are derived in §3. Numerical 
results are summarized in §4 and discussions about magnetic modes are summarized in §5. we conclude in §6. The details of 
the oscillation equations solved in this paper and suitable boundary conditions imposed at the stellar center and surface are 
given in Appendix. 


2 EQUILIBRIUM MODEL 

Although we are concerned with magnetized neutron stars, in this study, we consider the problem within the frame work of 
Newtonian ideal magnetohydrodynamics for the sake of simplicity. For the modal analysis of magnetized stars composed of 
the infinitely conductive fluid, in this paper, magnetic fields in equilibrium are treated as perturbations about non-magnetized 
and non-rotating star. As mentioned before, we ignore the deformation of the star due to the magnetic stress and use self- 
gravitating polytropic spheres for the matter distribution of the star. Then, the mass density and pressure are given by 

p = po0", J5 = po0"’^\ (1) 


where po and po are the density and pressure values at the center of the star, respectively, and n and 0 are the polytropic 
index and the Lane-Emden function, respectively. Imposed stationary axisymmetric magnetic fields are assumed to be the 
purely poloidal and dipole ones, given by 


Br = 2/(r) cos 6, Bg = — 




sin0, Ba, = 0, 


( 2 ) 


in which V • B = 0 is automatically satisfied. Here and henceforth, the spherical polar coordinate (r, 6, <j)) has been employed. 
The function /(r) in Eq. (2) is determined by the Ampere law, V x B = 47rji,^r sin^e,^, which leads to 


d}-2 (Jj, 


-47rjV, 


( 3 ) 


where is related to the toroidal current and needs to satisfy integrability conditions for the ideal magnetohydrodynamic 
(MHD) equations. In this study, we assume that = cop where co is a constant determined by the boundary condition at 
the surface of the star. Near the center of the star, the function / behaves as 

f = ao + 0{Y), (4) 


where tro is a constant determined by the boundary condition at the surface of the star. We assume that j^f, = 0 outside the 
star. Thus, the exterior solution is given by = pb/A, where pb is the magnetic dipole moment of the star. We 

determine the constants ao and cq so that the interior solutions / and d//dr are continuously matched with the exterior 
solutions and d/*'®’‘^/dr at the surface of the star. The magnetic field lines in the polytropes of index n — I and 1.5 are 
shown in Figure 1. We see that there exist field lines closed within the star. 


3 PERTURBATION EQUATIONS 

The governing equations for non-radial oscillations of magnetized stars are obtained by linearizing the ideal MHD equations. 
Since the equilibrium state is stationary and axisymmetric, the time and azimuthal dependence of the perturbed quantities is 
given by the factor exp[i((Tt -|- where m is the azimuthal wave number. The linearized basic equations that govern the 

adiabatic, non-radial oscillations of magnetized stars are written as 

- = -V4>' - iVp' + ^ ^ [(V X B') X B + (V X B) X B'] , (5) 

P' + V • (p^) = 0, (6) 


P 


p' 

Fip 



( 7 ) 



non-axisymmetric magnetic modes 3 




Figure 1. Magnetic field lines in the polytropic model of index n = 1 (left panel) and n = 1.5 (right panel). 


B' = V X X B) , (8) 

where (') indicates the Eulerian perturbation, and rA in equation 0 denotes the Schwarzschild discriminant defined as 
dlnp 1 dlnp 


rA — 


d In r El d In r ' 


( 9 ) 


and Fi = {d\np/dlnp)sA- For polytropes of the index n, the adiabatic exponent for the perturbations are assumed to be 
given by 


n + 1 


+ 7 


( 10 ) 


with 7 being a constant, for which rA = — 7 (dlnp/dlnr) = 'yV. The star may be called radiative for 7 < 0, isentropic for 
7 = 0, and convective for 7 > 0. For radiative stars, we have p-modes, whose oscillation frequency is proportional to ^— 7 . 
For simplicity, we employ the Cowling approximation, neglecting 

Because of the Lorentz force term in equation 0 , separation of variables for the perturbations is impossible between the 
radial coordinate r and the angular coordinate 9. We therefore expand the perturbations in terms of the spherical harmonic 
functions (j>) with different I’s for a given azimuthal index m. The displacement vector ^ is given by (see e.g., Lee 2005, 

2007) 


£rS7(r)y*7(0,^), 


( 11 ) 


J = 1 




j=i 




( 12 ) 


= E 

j = l L 


rHi (r)- 


1 d 


sm 9 o4> 1 a9 1 


and the vector B' is given by 

B:='|7%fo^.(r)y,7*(0,^), 

^ 3 3 


( 13 ) 


(14) 


B'e=Y. 


J = 1 




(15) 






(16) 


where Ij = |m| + 2{j — 1 ) and Ij = Zj■ + 1 for even modes, and Ij = |m| +2j — 1 and Ij = Ij — 1 for odd modes, respectively, 
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and j = 1, 2, 3,jimax- The Euler perturbations of the pressure and density are given by 

Jmax Jmax 

p' = PhY)yi^{s^4>), p' = (17) 

t=l 

For stars having a purely poloidal magnetic field, the oscillation modes are separated into even and odd modes. If the sets of 
the functions (^r, Bg) and {^g, B', B'A of a mode are respectively even and odd (odd and even) functions about the equator, 
we call the mode an even (odd) mode. In this paper, we usually use jmax = 12 to obtain solutions with sufficiently high-angular 
resolution. Substituting the expansions (11)-(17) into the linearized basic equations (O to ( 8 ), we obtain a finite set of coupled 
linear ordinary differential equations for the expansion coefficients Si . (r), T;/ (r), bf (r), bp (r), bT, (r), and rdbp (r)/dr, which 
we call the oscillation equations to be solved in the interior of magnetized stars. The oscillation equations obtained for the 
star magnetized with a purely poloidal field are given in Appendix A. The set of ordinary differential equations is solved as 
an eigenvalue problem of by applying boundary conditions at the center and surface of the star (see also Appendix A). 
Since the eigenvalue is a real number for the boundary conditions we use, the eigenvalue > 0 corresponds to stable and 
purely oscillatory modes of the frequency icr, and < 0 to unstable and monotonically growing modes of the growth rate 
rj = 

For a given jmax, we find numerous solutions to the oscillation equations. Most of them, however, are dependent on jmax. 
We have to look for solutions that are independent of jmax. 


4 NUMERICAL RESULTS 

We use polytropes of index n = 1 and 1.5 as a background model to calculate non-axisymmetric (m 7 ^ 0) oscillations of stars 
magnetized with a poloidal field. In this numerical study, magnetic modes, the eigenvalue cr^ of which is proportional to B|, 
are only the modes we discuss. We cannot correctly compute g-, /-, and p-modes of the magnetized stars (see § 6 ). We find 
both stable (oscillatory) magnetic modes with cr^ > 0 and unstable (monotonically growing) magnetic modes with a < 0. In 
fact, if we write cr = mi = ±pi with rj being a positive real number, the time dependence of the modes is given by exp(=p 77 t), 
and the modes with exp( 77 t) monotonically grow with time without bound where rj may be regarded as the growth rate. The 
modes with < 0 may correspond to magnetic instability. It is well known that the stars having purely poloidal magnetic 
fields are unstable and the energy of the field is dissipated quickly, that is, for several ten milliseconds (e.g., Markey & Tayler 
1973, van Assche, Goossens, Tayler 1982, Braithwaite 2007, Lasky et al. 2011; Ciolfi & Rezzolla 2012). 


4.1 Stable Magnetic Modes 

In Figure 2, we plot the eigenfrequencies a of stable magnetic modes that have no radial nodes of Si, for m = 1, 2, 3, and 4 
versus the Alfven frequency a a defined as 

a A = Bs/-J drrpc.BS, (18) 

where pc. is the central density, R is the radius of the star, and Bs = pbIR? is the magnetic field strength measured at 
the surface, and the frequencies a and a a are normalized by Hk = \/ GM / where M is the mass of the star and G is 
the gravitational constant. In this paper, we assume M = 1.4Mq and R = 10®cm, for which we have the ratio oaI^k ~ 
4.42 X 10“"^(Bs/10^®G). In this paper we use the central density pc to define the Alfven frequency cta. We may use instead 
the mean density p = M/( 47 rB^/ 3 ) for the definition. Since for polytropes p = —3pcCr^ = Cpc, where 0(C) is 

the Lane-Emden function and 0(Ci) = 0, we have C = Z/'Y ~ 0.304 and 0.167 for the indices n = 1 and 1.5. If we define 
(j*A = Bs/ A-KpR? , we have (j*a — (ja/'/C- Note also that Q,k = Y dnGp/Z oc \/Gp. We find that the eigenfrequency of 
the modes is isolated and proportional to the Alfven frequency cta, that is, Bs- This property confirms that the modes we 
obtained are discrete magnetic modes. Note that we find stable magnetic modes only for odd parity and cannot find stable 
magnetic modes of even parity. 

In Table 1, we tabulate the eigenfrequency a = o/Q-k of stable magnetic modes for the polytrope of index n = 1 and 
1.5 for Bs = 10^® G. Here, we have assumed that 7 = 0 . Note that the magnetic modes we can find for each value of m are 
those that have only a few radial nodes of the expansion coefficient Si, and it becomes difficult to find magnetic modes as 
the number of radial nodes increases. From Table 1, we find that for given m and number of radial nodes, the normalized 
oscillation frequency d of the magnetic mode of n = 1.5 is larger than that of n = 1. For a given m, the oscillation frequency 
gradually decreases as the number of radial nodes increases. This tendency of the oscillation frequency is similar to that found 
by Asai & Lee (2014) for axisymmetric (m = 0) toroidal modes of the stars magnetized with a purely poloidal magnetic field 
in general relativistic framework. In addition, for a given number of radial nodes, the larger the azimuthal wavenumber m, 
the larger the frequency of the magnetic modes. 

In Figure 3, we plot the eigenfunctions of an m = 1 magnetic mode that has no radial nodes of Si, for the n = 1 polytrope 
and Bs = 10^® G. We find that the eigenfunctions S, H, and T of this mode have large amplitudes in the the core. We also 
find that the horizontal H and toroidal T components show rapid spacial oscillations near the stellar surface, although the 
radial component S does not. This phenomena will be discussed in §5. Note that for isentropic (7 = 0) models, we can obtain 
stable magnetic modes even for magnetic fields as weak as Bs ~ 10^^ G. 
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Figure 2. Eigenfrequency a of the stable magnetic modes of odd parity for m = 1, 2, 3, and 4 versus the Alfven frequency cr^ for the 
n = 1 polytrope. 


Table 1. Normalized eigenfrequency a of the stable magnetic modes of odd parity for Bq = 10^^ G. 


n = 1 

m 


number of radial nodes 



0 

1 

2 

1 

0.007523 

0.007120 

0.006954 

2 

0.007943 

0.007464 

0.007216 

3 

0.008174 

0.007710 

0.007439 

4 

0.008310 

0.007887 

0.007618 

n = 1.5 

m 


number of radial nodes 



0 

1 

2 

1 

0.010304 

0.009975 

0.009817 

2 

0.010720 

0.010247 

0.010038 

3 

0.010961 

0.010485 

0.010229 

4 

0.011106 

0.010663 

0.010397 



Figure 3. Expansion coefficients xSi, xTi/ ^ and xHi as a function of a: = r/R for an m = 1 stable magnetic mode of odd parity for the 
polytrope with n = 1 for Bs = 10^^ G, where the solid lines, the long dashed lines, the short dashed lines, and the dotted lines are for 
the expansion coefficients associated with Ij (or from j = 1 to 4. The amplitude normalization is given by = 1 at the surface. 
Here, the frequency d = cr/f^K of the mode is 0.007523. 
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Figure 4. Spatial oscillation patterns ^rix^O) (left), 0) (middle), and ^^(r, 0) (right) of the m = 1 stable magnetic mode of Figure 
3, where the amplitudes are normalized such that max(|^j('r, 0)|) = 1. 


Table 2. Amplitude ratios between the components of the displacement vector of stable magnetic modes of odd parity for = 10^^ 


G 


u = 1 

m 

i«r"i /!«>""! 


|5fjmax|/||max| 

1 

1.109 

0.322 

9.833x10-® 

2 

0.787 

0.304 

1.236x10“® 

3 

1.234 

0.366 

1.408x10“® 

4 

1.121 

0.396 

1.595x10“® 

n = 1.5 

1 

2.552 

0.111 

6.954x10“® 

2 

0.897 

0.047 

3.358x10“® 

3 

1.124 

0.051 

5.376x10“® 

4 

1.365 

0.103 

7.191x10“® 


Assuming 0 = 0, we may define the spacial oscillation pattern $ of the displacement vector $ as 

iAr,0) = = ( 19 ) 

where j = r, 6, (j), and the patterns ij{r,9) for an m = 1 stable magnetic mode are shown in Figure 4, where the vertical 
2 -axis given by x = 0 is the symmetry axis, and the amplitudes are normalized such that max(|^j (r, 0)|) = 1 for j = r, 6, (j). 
Since the mode is an odd mode, as shown by the figure, the patterns LYj S'lid 9) are antisymmetric about the equator 
given by 2 = 0, while ^e(r, 0) is symmetric. The oscillation patterns have large amplitudes along the symmetry axis. The ij) 
component of ^ shows a pattern that reflects the existence of magnetic fields closed within the interior. 

In Table 2, we tabulate the ratios S'lid for the nodeless magnetic modes of the n = 1 and 

n = 1.5 polytropes, where = max(|^j(r,0)|), and we have assumed 7 = 0 and Bq = 10^® G. From Table 2, we find 

that the relative amplitudes and are rather insensitive to the azimuthal wavenumber m, and that 

the components L and dominate The pressure perturbation SU = p'{r,9, 4> = Q)/pgr has almost negligible amplitudes 
compared to the displacement vector. The oscillation patterns of the nodeless mode are almost the same for the n = 1 and 
n = 1.5 polytropes although the ratios and for n = 1 are slightly larger than those for n = 1.5 as 

indicated by Table 2. For the magnetic modes that have non-zero radial nodes of Si,, the tendencies of the amplitude ratios 
are found similar to those for the nodeless magnetic modes, and the oscillation patterns look quite similar for the polytropes 
of n = 1 and 1.5 

As \m\ increases, the oscillation amplitudes tend to be confined in the envelope region away from the symmetry axis, as 
exemplified by the spatial oscillation patterns of an m = 5 unstable magnetic mode shown in Figure 8. 
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Table 3. Growth rate r} = of the monotonically growing magnetic modes of even and odd parities for Bg = 10^® G 


n = 1 




even parity 




odd parity 


m 


number of radial nodes 


m 


number of radial nodes 



1 

2 

3 


1 

2 

3 

1 

0.000961 

0.000580 

0.000405 

1 




2 

0.001557 

0.000977 

0.000709 

2 

0.003562 



3 

0.001964 

0.001273 

0.000949 

3 

0.004435 

0.002038 


4 

0.002260 

0.001505 

0.001144 

4 

0.004879 

0.003088 


5 

0.002488 

0.001695 

0.001452 

5 

0.005147 

0.003685 

0.002008 

6 

0.002676 

0.001859 

0.001452 

6 

0.005326 

0.004086 

0.002728 

7 

0.002852 

0.002110 

0.001889 

7 

0.005454 

0.004377 

0.003209 

8 

0.003096 

0.002616 

0.002069 

8 

0.005548 

0.004598 

0.003566 

n = 1.5 



even parity 




odd parity 


m 


number of radial nodes 


m 


number of radial nodes 



1 

2 

3 


1 

2 

3 

1 

0.000898 

0.000430 


1 




2 

0.001601 

0.000804 


2 




3 

0.002114 

0.001121 


3 




4 

0.002482 

0.001386 


4 

0.006170 

0.003683 


5 

0.002769 

0.001618 


5 

0.006510 

0.004506 

0.002007 

6 

0.003328 

0.002611 


6 

0.006734 

0.005047 

0.003147 

7 

0.004089 

0.002895 


7 

0.006892 

0.005435 

0.003849 

8 

0.004659 

0.003053 


8 

0.007009 

0.005727 

0.004354 


4.2 Unstable Magnetic Modes 

Unstable magnetic modes are found both for even parity and for odd parity. If we define the growth rate ty > 0 such that 
cr = ±? 7 i, rj is almost exactly proportional to the field strength Bg, particularly for the case of isentropic (7 = 0) models. In 
Table 3, we tabulate the normalized growth rate fj = tj/Qk of unstable magnetic modes for the polytropes of index n = 1 and 
1.5, where we have assumed 7 = 0 and Bg = 10*^® G. Note that unstable magnetic modes are found even for magnetic fields 
as weak as Bg ~ 10*^^ G. From Table 3, we find that the growth rates fj of the unstable magnetic modes for the n = 1 and 
n = 1.5 polytropes are quite similar, and that it becomes difficult to obtain unstable magnetic modes as the number of radial 
nodes increases, particularly for n — 1.5. Between m — 1 and m = 8, we find unstable magnetic modes of odd parity only for 
m ^ 2 for n — 1 and for m ^ 4 for n = 1.5. For a given m, the growth rate fj gradually decreases as the number of radial 
nodes of Sq increases. On the other hand, for a given number of radial nodes, the larger the azimuthal index m, the larger 
the growth rate of the unstable modes. Using the normalized growth rate fj, the growth time scale Tg = l/rj may be given by 

Tg = 7.33 X 10~®/?7 sec, (20) 

where we have assumed M = \AMq and R = 10®cm. For a typical value of 7 ~ 10“®, the growth time scale may be 
Tg ~ 5 X 10“^sec, which is consistent with the results, e.g., by Lasky et al. (2011) and Giolfi & Rezzolla (2012). 

Figure 5 shows the eigenfunctions of an m = 1 unstable magnetic mode of even parity that has one radial node of Sq for 
the n = 1 polytrope for Bg = 10^® G, where the amplitudes are normalized by Tii^ = 1 at the surface. We find that 5';^ has 
large amplitudes in the core, while in the envelope. The function 7/q has large amplitudes both in the core and in the 
envelope of the star. Note that the first components of S, T, and H are dominating the other components. 

The oscillation patterns ^r{r,9), ^e(r,9), and ^4,{r,9) for the m = 1 unstable magnetic mode are given in Figure 6 . Since 
this mode has even parity, and are symmetric about the equator, while is antisymmetric. The patterns has large 
amplitudes along the equator, but the amplitudes of the 9 and 0 components of ^ are confined to the polar regions. The 
existence of closed magnetic fields is recognized in the oscillation pattern of ^r- 

In Table 4, the ratios S'lid tabulated for the unstable magnetic modes of the n = 1 and 

n = 1.5 polytropes for Bg = 10^® G. From Table 4, we find that relative amplitudes decrease 

with the azimuthal wavenumber m, and that dominates the other components. The pressure perturbation SU have negligible 
amplitudes compared to the displacement vector as in the case of stable magnetic modes. The oscillation patterns are almost 
the same for the n = 1 and n — 1.5 polytropes. 

Figures 7 and 8 show the eigenfunctions and spatial oscillation patterns of an m = 5 unstable magnetic mode of odd parity 
that has one radial node of Sq for the n — 1 polytrope for Bg = 10^® G. The toroidal components Ti^ of the displacement 
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Figure 5. Same as Figure 3 but for an m = 1 unstable magnetic mode of even parity with the growth rate 77 = 0.00096. 



Figure 6. Same as Figure 4 but for the m = 1 unstable magnetic mode of Figure 5. 


vector ^ look similar to those of the m = 1 unstable mode of even parity in Figure 5, but the amplitudes of Si^ and Hi. are 
confined in the envelope region, in contast to those of the m = 1 magnetic mode. This amplitude confinement to the envelope 
region is clearly seen in the spatial oscillation patterns shown in Figure 8. It is also interesting to note that the mode has 
negligible amplitudes along the symmetry axis. 


Table 4. Amplitude ratios between the components of the displacement vector of the unstable magnetic modes of even parity for 
Bs = IQl® G 


n = 1 

m 


||max|/||max| 

|5fjmax|/||max| 

1 

17.05 

319.1 

3.462x10"'^ 

2 

10.95 

120.5 

3.706x10-“* 

3 

9.183 

72.64 

3.137x10-“* 

4 

7.996 

49.24 

2.538x10-“* 

n = 1.5 

1 

31.98 

685.1 

1.471x10-“* 

2 

14.39 

159.9 

1.166x10-“* 

3 

7.662 

63.72 

7.266x10-5 

4 

10.66 

69.91 

1.082x10-“* 
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Figure 7. Same as Figure 3 but for an m = 5 oscillatory magnetic modes of odd parity with the growth rate rj = 0.005147. 



0 0.5 1 


x/R 



Figure 8 . Same as Figure 4 but for the m = 5 unstable 


magnetic mode of Figure 7. 


5 DISCUSSION 


Magnetic modes depend on 7 as shown by Figure 9, in which the frequency ct of a stable m = 2 magnetic mode of odd 
parity with no radial nodes of Si^ and the growth rate rj/aA of an unstable m = 2 magnetic mode of even parity are plotted 
as a function of for the n = 1 polytrope and for Bs = 10^® G. As I 7 I increases, the frequency cr gradually increases 

(decreases) for 7 < 0 (for 7 > 0), and we find no magnetic mode for \/py[/'^A ^ 10. On the other hand, the growth rate rj 
decreases (increases) as I 7 I increases for 7 < 0 (7 > 0), and we hud no unstable solutions beyond > 4. For radiative 

stars with 7 < 0, the buoyant force tends to stabilize the magnetic instability. It is important to note that when we normalize 
the eigenvalue a (or r\) and the Brunt-Vaisala frequency N oc \/py[ in terms of cta, the relation between aja a (or rjjaA) and 
| 7 y^/(TA in Figure 9 is almost independent of the magnetic held strength Bs- 

Using a dispersion relation derived by Lee (2010) for the oscillation of magnetized stars, we try to explain the rapid 
spatial oscillations of the expansion coefficients H and T in the surface layers for stable magnetic modes. For a non-rotating 
and isentropic star, the dispersion relation may be given by 


— pq cos a {Rk) + q cos a \ 2- + 1 ] a (Rk) + q 


q 2 2 
p cos a — 


<7 

+a‘^ \a* — + q (Rkn)^ 0^ sin^ -!/>] = 0 , 


-I- 1 -b cos a a 


a {Rkf 


where 

aP _ B^! Ait p ^2 

{RQ.kY' ^ 


2 / 2 
9 /a 
02 
“/c 


p GMr k ■ B . 

Fi-, g= ——, cosa=—, smV’ 
p r2 kB 


(fc X B), 

knBH 


( 21 ) 


( 22 ) 


and k = |fe| and B = |B|. When we employ local cartesian coordinates {x,y,z) and assume that the 2 -, x-, and ^-directions 
at a point in the interior are respectively along the r-, 9-, and ^-directions at that point, we have B^ = 2/cos0, By = 0, 
Bx = — (2/-b d//dlnr) sin0, kn = yfcf + fcj, Bh = \Bx\, and sinTp = —kyBx/kH\Bx\- Assuming a polytrope of index n = 1 
with the mass M = IAMq and the radius R — 10®cm and the held strength Bs = 10^®G, we can compute the quantities p, q, 
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Figure 9. Eigenfrequency crjaA (left) and the growth rate rijaA (right) of the m = 2 magnetic modes versus the for the 

n = 1 polytrope. 


\/R 



r/R 


r/R 


Figure 10. Local wavelengths X/R = 2iT/{Rkz) derived using the dispersion relation El for the n = 1 polytrpoe of M = IAMq and 
R = 10®cm and for Bg = 10^®G, where the black lines and red lines are for cos0 = 0.5 and 0, respectively, and the solid lines, dashed 
lines, and dotted lines are for cosq = 0.5, 0.1, and 0.01, respectively. Here, we assume = 10“® (left) and = 10““^ (right). 


, a^, y as a function of r. If we further assume for simplicity kxR = kyR = 1 for \m\ ~ ~ 1 in the surface regions, we may 
solve equation dni) for (Rk)^ to obtain (Rkz)^ for given values of a^, cos 9, and cos a. Figure 10 shows the local wavelengths 
X/R — 2Tv/{Rkz) computed using the largest positive solution (Rkz)^ in the outer envelope for several sets of parameters 
((j^, cos 6^, cos Of). As shown by the figure, the waves propagating across the held lines, which may be provided by the closed 
held lines near the surface, tend to have short wavelengths, which become even shorter toward the surface. The wavelengths 
of the rapid spatial oscillations of the functions H and T near the surface become comparable to those expected for the case 
of cosa ~ 0.01. For low frequency modes we may expect fe • ^ ~ 0, and hence <C 1 for \kz\/kH 1- In this case, 

rapid spatial oscillations in the radial components S in the surface layers may not be conspicuous. 


6 CONCLUSION 

In this paper, we have calculated non-axisymmetric (m yf 0) oscillations of neutron stars magnetized with purely poloidal 
magnetic helds, where we have used polytropes of index n = 1 and 1.5 as background neutron star models. We have found 
stable (oscillatory) magnetic modes (cr^ > 0) of odd parity and unstable (monotonically growing) magnetic modes {a^ < 0) 
of both even and odd parity. The frequency a of the stable magnetic modes and the growth rate r] = \/—a'^ of the unstable 
magnetic modes are proportional to the magnetic field strength Bs = fJ-b/R^ measured at the surface, if the effects of buoyancy 
in the interior is negligible. For a given m, the frequency a and the growth rate rj decrease as the number of radial nodes 
of the eigenfunctions increases, which may indicate an anti-Strumian property of the problem. We have found that the non- 
axisymmetric magnetic modes are affected by stratihcation in the interior of the star, which is parametrized by using a 
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parameter 7 in this paper. For radiative stars (7 < 0), the eigenfrequency a of the oscillatory magnetic modes gradually 
increases as I7I increases, while it decreases with I7I for convective stars (7 > 0). For the unstable magnetic modes, on the 
other hand, we have found that stable stratihcation with 7 < 0 reduces the growth rate rj of the magnetic modes, while the 
convectively unstable stratification enhances the growth rate as I7I is increased. It is also found that the growth rate rj tends 
to increase as the azimuthal number increases. We note that Lasky et al. (2011) obtained strong instability of magnetized 
star having purely poloidal magnetic field, especially for m = 4. It has been analytically shown that purely poloidal magnetic 
helds having closed held lines inside the star are unstable in the limit of m ^ 00 (see, e.g., van Assche, Goossens, & Tayler 
1982, Markey & Tayler 1973). 

As mentioned in §4, we fail to obtain /- and p-modes of the magnetized star by using the present numerical code. At hrst 
glance, this fact might imply that the present numerical code and/or formulation have some problem. However, the situation 
at hand is not very simple, as we will discuss below. Let us assume that in the limit of \B\ —>■ 0, there exist eigensolutions for 
magnetized stars given by 

(7 = no + 0(di), ^ = ^0 + O(di), p'=Po + 0(di), p'= p(, + O(CTi) , (23) 

where ao is an eigenfrequency and ^g, pg, and pg are its eigenfunctions for non-magnetized stars. Clearly, Eq. (23) is an 
eigensolution of Eqs. (5)-(7) in conjunction with the surface boundary condition Ap(r = i?) = 0 in the limit of |B| —>■ 0. 
From Eq. (8), we obtain 

J = vx(e„x.^)+o(7). (24) 

b' b' 

However, in general, this does not satisfy the surface boundary condition that guarantees —>■ 0 as r —>■ 00. Thus, 

Bs Bs 

Eqs. (23) and (24) cannot be an eigensolution for magnetized stars in the limit of \B\ —>■ 0. In other words, oscillation modes 
that exist in non-magnetized stars need not be present in weakly magnetized stars. From physical point of view, however, 
it is plausible to expect that /- and p-modes exist in magnetized stars. This apparent contradiction may be resolved if we 
admit that in the normal mode analysis such as we carry out in this paper the Lorentz force term in Eq. (5) must not vanish 
even in the limit of |B| —>■ 0. This occurs if V x B' ex is satisfied somewhere inside the star when |B| —>■ 0. That is to 
say, B' needs to show rapid spatial oscillations somewhere inside the star to keep the Lorentz terms finite as \B\ —>■ 0. If 
this is the case, /- and p-modes of magnetized stars have to have slightly different eigenfrequencies and eigenfunctions form 
those of non-magnetized stars even in the limit of \B\ —>■ 0. In the present numerical code, it will be difficult to treat the 
terms related to V x B' oc accurately because we use a standard second-order accuracy Hnite-difference method for the 
radial direction and a spectral method for the 0-direction, and some special technique will be required to evaluate the terms 
related to V x B' oc d)( properly. Actually, we find that we can obtain ‘/’- and ‘p’-modes if we largely reduce the number 
of mesh points distributed in the interior for integration. This reduction in the number of mesh points may be equivalent to 
averaging the rapid spatial oscillations of B' over a length scale much larger than the wavelengths of the spatial oscillations, 
which avoids numerical difficulty associated with rapid spatial oscillations of B' . However, we cannot consider that ‘/’- and 
‘p’-modes thus obtained are correctly computed normal modes. 

It is also important to note that, for given values of m and Bs, we hnd only one magnetic mode sequence along the 
number of radial nodes of the eigenfunction. This result is different from that for the axisymmetric toroidal magnetic modes of 
stars with a poloidal field as discussed by Lee (2008) and Asai & Lee (2014), who obtained several mode sequences, differing in 
the surface oscillation patterns, for given m and Bs- If discrete magnetic modes can exist only in the gaps between continuum 
frequency spectra, the difference in the distribution of continuum frequency spectra between axisymmetric (m = 0) toroidal 
modes and non-axisymmetric (m 7^ 0) spheroidal modes, for given m and Bs, might lead to the difference in the distribution 
of discrete magnetic modes between the two cases. 

The present analysis is a part of our study of the oscillation of magnetized stars. It is well known that a purely poloidal 
magnetic field conhguration is unstable as exemplified in this paper. Therefore, it may be difficult to observe long-lived 
magnetic oscillations as QPOs. However a mixed poloidal and toroidal magnetic held conhguration such as twisted-torus 
magnetic held (e.g., Braithwaite & Spruit 2004; Yoshida & Eriguchi 2006; Yoshida, Yoshida & Eriguchi 2006; Ciolh et al. 
2009) could be stable. Thus, it is interesting to investigate stability and oscillation spectrum of the star having such a mixed 
magnetic held conhguration. In the presence of both poloidal and toroidal held components, toroidal and spheroidal modes 
are coupled even for axisymmetric perturbations, which inevitably makes the analysis difficult. 

As important physical properties inherent to cold neutron stars, we need to consider the ehects of solid crust and those 
of superhuidity and super-conductivity of neutrons and protons on the oscillation modes. It is believed that neutrons become 
a superhuid both in the inner crust and in the huid core while protons can be superconducting in the core. For example, 
if the huid core is a type I superconductor, magnetic helds will be expelled from the core region bacause of the Meissner 
ehect, and hence conhned to the solid crust (e.g., Colaiuda et al. 2008; Sotani et al. 2008). However, a recent analysis of the 
spectrum of timing noise for SGR 1806-20 and SGR 1900-1-14 has suggested that the core region is a type II superconductor 
(Arras, Gumming & Thompson 2004). If this is the case, the huid core can be threaded by a magnetic held and hence the 
frequency spectra of oscillation modes will be affected by the superconductivity in the core (e.g., Colaiuda et al. 2008; Sotani 
et al. 2008). To investigate the oscillation of magnetized stars as normal modes, taking account of the effect of superhuidity 
or superconductivity, will be one of our future studies (see, e.g., Glampedakis, Andersson & Samuelsson 2011; Passamonti & 
Lander 2013, 2014; and Gabler et al. 2013b). 
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APPENDIX A: PULSATION EQUATIONS FOR THE MAGNETIZED STAR WITH PURELY 
POLOIDAL MAGNETIC FIELDS 


To describe the master equations concisely, it is useful to introduce the following column vectors composed of the expansion 
coefficients for the perturbation quantities: the vectors S, H, T, b®, 6^, b^, and 5U, defined by 


{S)i=Si^, {Th=T, {b%=bP,, {b’^L^bT., (SU), = 

J 3 j ^ pgr 


(Al) 


where (X)j denotes the j-th component of the vector X and g — GM{r)/Y is the gravitational acceleration. Here, Ij = 
|m| + 2{j — 1) and Ij = Ij + 1 for even modes, and Ij = |m| + 2j — 1 and Z' = Ij — 1 for odd modes, respectively, and 
j — 1, 2,3,..., jmax. The perturbed continuity equation (5) and the perturbed Euler equation (4) are then reduced to 


dS 


r 


dr 



S-^SU + l\oH, 

i 1 


(A2) 


_ + rA)S + ^V{1 -rA- U)SU 

dr rj rj 


+Co 


(d/ + 2)r^ - (d/ + 2)b® + (dV + 6d/ + 4)b" 


+ m 


(d/ + 2)r^|^ + {(Y f + 6d/ + 4)b^ 


= 0 , 


(A3) 


Lsu = ^Cid2(d/ + 2)Ar'BiS - + Aq + j^cid'^Ar' 


+ 


rf 

27rr ip 


(dV + 4d/)Ao ^Bo - Ao'Ao + f + 4d/ + 2)Ao ^Ar^ | 


+^^(Ao ^Ao - m^Ao ^Ar^)b" + -r^^m 


ttFip 


d/Ao'Ar'BiAo + Ao'Ar'Ci 


27rrip 

b^ + 


i(d/ + 2)l + 2Ao' 


rf 

27rr ip 


Aq AoLor 


_ 

dr 


(A4) 


Air^^;^ = —^^Uckt^AiT + ^m(d^/ + 4d/ + 2)b® — 2mb^ + ("d/BiAo + Ci) b^ — . 

dr rj 2 \ / dr 

Substituting 5U given in (A4) into (A3), we may obtain 

Lor^ (r^) = (ci<t" + rA) Aq ^Ao5 + ^mUcid^Ao'Ar'AiT + | - i (d/ + 2) Aq ^Cq - (dp - rA + 4) I 

+ [-d®/ + d^/ (dp - rA - 4) + 4d/ (dp - rA + 1)] Aq ^Bq - ^ [d®/ - d^/(dp - 4) 

— 2d/(2dp +1) — 2(dp + 2)] m Aq A^^ A [(d/ + 1)Aq Cq + Ai + 2 Aq BqAi + 2(dp — rA)l 


+ 


(A5) 


--(d7 + 4d/ + 2)m"Ao'ArWi - (d7 + 2df + 2dp)m"Ao''Ar 
I - m|2(d/ - dp + rA + l)Ao-' + i[-d/(dp + 2) + (d/ + 2)rA - 2(dp + 1 )]Ao-1Ao| 
-m [(d7 - dpdf + 2d/)Ao-7r^BiAo + (d/ - dp + l)Ao 7r7i] |b^ 


+ 


-(d/ + 2)Ao 7o + (dp - rA - 1)1 + m"(d/ - dp + 3)Ao 


r^ - ^(2 + rA)Uci7Ao ^AoH 
dr rj 


27vp^^ d// 27ip -2.-1--1. dT 

-Aq Aor—-^-j-mFciu Aq Ai Air—-m 

rj dr rj dr 


i(d/ + 2)Ao7r7i + (d/ + rA + 2)Aoi r^. (A6) 


The perturbed induction equation (7) and the perturbed Gauss’s law for magnetic fields are reduced to 
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Lor— = -m dV + 2 + — d/ - 2 1 - — Aq'A r'Bi + -Aq ^ S + (d/ + 2) m^Ap ^Ar' + ^Aq'C o T 


mAo-'Ar'Aib" + — Ao ^Aob^ - m (d/ + 2) (^Aq ^ + |Ao "Af^CiJ H 

+m (df + 2 ) ('ao'A r'Bi + iAo ^SU, (A 8 ) 


k—1 I ■*'A—1 a—1/ 


Lir— = dV + 2 + — d/ - 2 1 - — Ar'Bi + -m^Ar^AQ-' S - m (d/ + 2) A^' + -Ar'Ao-'Co T 


+—Ar'Aib" - —mAr'Ao'Aob^ + {df + 2) (^m^Ai-^Ao ^ j H 

- {df + 2) f A^^Bi + im^Ar^Ao ^SU, (A9) 


r^ = -3b® + Aib". 
dr 


dlnp M / r 

^ "" “dlnr ’ ^ M{r) \r) ’ ^ 


dlnp dW , 2 . 

d In r ’ d In r ’ / dr^ ’ / dr® ’ 


and d = ol{GMIB?fl‘^ is the frequency in the unit of the Kepler frequency of the star. The quantities Ao, Ai, Bo, Bi, Co, 
Cl, Ao, Ai, Lo, and Li denote the matrices defined as follows: 

Ao = Co + QqAi, Ai = Ci + QiAo, Bo = Qo + 2 ^'^’ d" 2 ^^’ 

Co = Co(Ai — 21), Cl = Ci(Ao — 21), Ao = Ao — 21, Ai = Ai — 21, 

Lo = I - m^Ao’^Ar\ Li = I - m^Af^Ao \ (A12) 

The matrices Qo, Qi, Co, Cl, Ao, and Ai are defined as follows: 

For even modes, 

(Qo)ii = J/7+1, (Qo)i+i,i = J7+2, (Qi)ii = J7+i, (Qi)i,i+1 = J7+2, 

(Co)ii = -{h + 2)jr+i, (Co)i+i,i = {h +1)jr+ 2 , (Ci)ii = Jz7+1, (Ci )„+1 = -{I, + 3)jr+ 2 , 


(^o)ii — lj{lj + 1), (Ai)jj — {Ij + l){lj + 2), 
where Ij = |m| + 2j — 2 for j = 1, 2, 3,...., jimax, and 
im _ r {lj+m){lj-m)Y^^ 


L(2/i-l)(2Zi + l)J ■ 

For odd modes, 

(Qo)ii = J7+1, (Qo)„+i = J7+2, (Qi)ii = J7+1, (Qi)i+i.i = J7+2, 

{^o)jj = (Co)j,j+i = —{Ij + 3) J 7 + 2 , (Ci)ii = —{Ij + 2) J(7+i, (Ci)j+ij’ = {Ij + l)Ji7+2, 


{^o)jj — {Ij + l)(^j + 2), (Ai)jj — lj{lj + 1), (A15) 

where Ij = \m\ +2j -1 for j = 1, 2, 3,...., jmax- 

From the equations given before, we see that non-axisymmetric pulsations of the magnetized star with purely poloidal 
magnetic fields may be described by a system of 6jmax-th order ordinary differential equations. In this study, we chose the 

o Tj '~p d 

column vectors S, T, b , b , b , and r — b as dependent variables. If the dimensionless vector variables, defined by 

dr 

Vi = S, 1/2 = T, ^3 = h^, ^4 = a", yg = h'^, (^^16) 

where A* = [i?//(0)]b* for i = S,H,T, are introduced, the master equations for stellar pulsations are schematically written 
by the coupled first-order differential equations, given by 

r^=:Fiiyi-hAoiF-^<5I7 
dr i 1 

= (-^11 — ^11 + Aof2l)yi + ( —f 12 + Aof22)2/2 + { — £13 + Ao£23)y3 — f 142/4 “ ^152/5 “ ^162/6i (^17) 
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= ^212/1 + ^222/2 + -^242/4 + -^252/5 + ^21^SU + Q22H 
= {^21 + ^21^11 + ^22^21)2/1 + (^22 + ^21^12 + ^22^22)2/2 “I” (^21^13 + Q 22 ^ 2 z ) y ^ 

+ (^24 + ^21^14)2/4 + (-^25 + ^21^15)2/5 + ^21^162/65 


^ d2/3 

dr 


-82/3 + A12/4, 


/?/4 

dr 


= ^6- 


dr 


^522/2 + -^532/3 + -^542/4 + -^552/5 + -^562/61 


r — ^6lPl + -^622/2 + -^632/3 + -^642/4 + -^652/5 + -^662/6 + ^62-ff + ^63’'+ Qq^T-^ 
= [-^61 + ^62^21 + ^ 63(^31 + £ 37^21 + fsS^^ll) + 064(-?^21 + 021^^11 + ^22^21)] 2/l 

+ [-^62 + Qm^ll + ^63(^32 + f 37^22 + 38^12) + S64(-^22 + ^21^12 + 022^22) + ^65-^52] 2/2 

+ [^63 + 062^23 + ^63(^37^23 + ^38^13) + ^64(^21^13 + ^22^^23) + ^65-^53] 2/3 
+ [-^64 + ^63(^34 + f 38^14) + ^64(-^24 + ^21^14) + ^65-^54] 2/4 

+ [-^65 + ^63(^35 + ^38^15) + ^64(-^25 + ^21^^15) + ^65-^55] 2/5 

+ (-^66 + ^63^38^^16 + 064^21^16 + ^65-^56) 2/6i 


and the algebraic relations, given by 

y 

—( 5/7 = £iiyi + f i2y2 + 132/3 + ^142/4 + ^152/5 + ^^i 62 / 6 > 


H = £ 21^1 + f 22y2 + £23y3, 

= £ 312/1 + £ 322/2 + £ 342/4 + £ 352/5 + £3777 + £38^(5/7 

= (£31 + £37£21 + £38£ii)2/i + (£32 + S-37^22 + £38£l2)2/2 + (£37£23 + £38£l3)2/3 

+ (£34 + £38£i4)2/4 + (£35 + £38£i5)2/5 + SssSwU^- 
The coefficient matrices appearing in (A17)-(A25) are defined by 
£11 = Vccid" (d/ + 2) A^^Bi, £12 = -mVacia^ (Ar^ + Aq ^Ar^Ai) , 


£13 = —.Vccia^fi^A + il-VccxYi 


2 / 


2px 


(dV + 4 d/) Ao'Bo - Ao'Ao + -rrd {<Yf + 4 d/ + 2) Aq'A r' 


(A18) 

(A19) 

(A20) 

(A21) 


(A22) 

(A23) 

(A24) 

(A25) 


£i 4 = -^VgciujI (Aq-'Ao - 
px 


f -2 

£15 = -—VaciLOAm 
2px 


i (d/ + 2 ) I + 2 Ao ' + d/Ao'A j^'BiAo + Aq'Aj:'Ci 


£i6 = -^VgCIOIaAo 'AoLo, 
2px 


£21 = (d/ + 2 )Ar'Bi, £22 = -mAr\ £23 = ^Af', 

2/ 

£31 = [dV + (2 + 1 /g) d/ - 2 (1 - Vg)] L7' (^A^'Bi + im^Ar'Ao , 

£32 = -m (d/ + 2 ) L 7 ' (^Ar' + iA^'Ao'Co) , £34 = ^Lr'A^'Ai, £35 = -^mL^'Ar'Ao'A o, 

\ / 2 f 2/ 

m^Ar'Ao , (A26) 

7^11 = (T7g - 3) I, 7^21 = -m [dV + (2 + Vg) d/ - 2 (1 - Vg)] Lq ' (^Aq 'Ar' Bi + Iaq, 


£37 = (d/ + 2 ) Lr^ 


^Ar'Ao-' + ^Ar'ci^, £38 = - (d/ + 2 ) Lr' (^Ar'Bi + i 
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= (d/ + 2) m"Ao + ^Aq "Co , J^24 = -.mLo "Ao'"Ar"Ai, J^25 = ^Lo'"Ao'"Ao, 

\ ^ / 2/ 2/ 

:F52 = -2^^Ar"Ai, :F53 = i (dV + 4d/+ 2) mAr", =-2mA^\ jFgs = A^" fd/BiAo + Ci) , 

f lUa 2 V / 

^56 = —m/\i ", = 2 


-r icid" + rA) Lo "Ao'"Ao, JFgz = 4^^mLo-'"Ao "Ai-"Ai, 

/ 


:^63 = Lo^ (d/ + 2) Ao^"Co - (dp + 4 - rA) I + [-dV + d"/ (dp - 4 - rA) + 4d/ (dp + 1 - rA)] Ao^"Bo 

-i [dV - dV (dp - 4) - 2d/ (1 + 2dp) - 2 (2 + dp)] m^Ao "Ar"|, 


= Lo " 


(d/ + 1) Ao "Co + Ai + 2 Ao'"BoAi + 2 (dp - rA) I - - (dV + 4d/ + 2) m^Ao "Ar"Ai 


- (dV + 2d/ + 2dp) m^Ao "AjT" 


:^65 = Lo " I - m 12 (d/ - dp + 1 + rA) Ao " + i [-d/ (dp + 2) + (d/ + 2) rA - 2 (dp + 1)] Aq " Ao | 

-m [(dV - dpd/ + 2d/) Ao "Ar"BiAo + (d/ - dp + 1) Aq "Ar"Ci] |, 


•^66 — Lo 


I (d/ + 2) Ao "Co + (dp - 1 - rA) I + (d/ - dp + 3) m^Ao "A^" 


(A27) 


^21 = m (d/ + 2) Lo " (^Ao "Ar"Bi + Iaq "j , ^22 = -m (d/ + 2) Lq " (^Aq " + ^Aq "Ar"Ci ) , 

^ -2 -'—2 -'—2 

062 = -2^^ (2 + rA) Lo-"Ao "Ao, 063 = -2^^Lo "Aq "Aq, 064 = 2^^mLo-"Ao "Ar"Ai, 

f A f f ^A 


(d/ + 2) Ao "Ar"Ci + (d/ + 2 + rA) Aq " 


^'es = -niLo 
where 

V - — ^ - L f - /ill A - ^ ,2 _ /^(O) p _ /^(o) 

Ti’ W ^ /(O)’ ^ p(0)’ ^ 7rp(0)i?2’ P 47rp • 

Using (A17)-(A25), we can summarize our master equations as follows: 


dr 


( Vi \ 


( Vi \ 

V2 


V2 

Vs 

= A 

Vs 

Vi 


V4 

Vs 


Vs 

\ Vs / 

\ Vs J 


(A28) 


(A29) 


(A30) 


where ,4. is a matrix derived through complicated calculations. We assume that no electric current flows outside the star. 
Therefore, the magnetic field perturbation regular at inhnity needs to satisfy the relations, given by 


= 0 , = 0 , for r > R , 

where (L^)ij = (Z) + l)Sij. The surface boundary conditions assumed in this study are given by 
[[AB]] = 0, Ap=0, 


(A31) 

(A32) 


where for the physical quantity Q, AQ denotes the Lagrangian change in Q and [[Qj] = lim{Q(7? — e) — <5(77+ e)} for e > 0. 

€—>■0 


Note that the mechanical surface boundary condition is usually given by 


a|. + T|B|“ 


term in the left-hand side of this equation automatically vanishes due to the conditions of [[AB] 
assumed in this study. For the magnetized star models employed in this study, the electric current vanishes at the surface of 


= 0. However, the second 
= 0 and [[Bj] = 0, which are 
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the star. We then have [[AB]] = [[B']] = 0, from which Eq. (IA31I) has to be satisfied at the surface of the star. The condition 
Ap = 0 is explicitly written by 

6U-y^=0. (A33) 

The boundary conditions at the stellar center are the regularity conditions for the eigenfunctions y^- y^. We adopt a nor¬ 
malization condition Tir^ {R) = 1 at the stellar surface. 
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